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Abstract 

The time-scale over which and modality by which young stellar objects (YSOs) 
disperse their circumstellar discs dramatically influences the eventual formation and 
evolution of planetary systems. By means of extensive radiative transfer (RT) mod- 
elling, we have developed a new set of diagnostic diagrams in the infrared colour-colour 
plane (K-[24] vs. K-[8]), to aid with the classification of the evolutionary stage of 
YSOs from photometric observations. Our diagrams allow the differentiation of sources 
with un-evolved (primordial) discs from those evolving according to different clearing 
scenarios (e. g. homologous depletion vs. inside-out dispersal), as well as from sources 
that have already lost their disc. Classification of over 1500 sources in 15 nearby star- 
forming regions reveals that approximately 39 % of the sources lie in the primordial disc 
region, whereas between 31 % and 32 % disperse from the inside-out and up to 22 % of 
the sources have already lost their disc. Less than 2 % of the objects in our sample lie in 
the homogeneous draining regime. Time-scales for the transition phase are estimated 
to be typically a few 10 5 years independent of stellar mass. Therefore, regardless of 
spectral type, we conclude that currently available infrared photometric surveys point 
to fast (of order 10 % of the global disc lifetime) inside-out clearing as the preferred 
mode of disc dispersal. 
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1 Introduction 



The lifetime and modality for the dispersal of protoplanetary discs around newly formed 
low mass stars (approximately solar mass or lower) is a key parameter that influences the 
formation and evolution of eventual planetary systems. The classical picture that emerged 
from the last decade of photometric observations, mainly carried out with the Spitzer Space 
telescope is that of disc evolution being described by two different timescales. The first 
timescale could be defined as a 'global timescale', i. e. the total time it takes from a star 
to go from disc-bearing to disc-less, and a 'dispersal timescale', i.e. the time it takes for a 
disc to disappear once dispersal ha s set in. Glo b al tim escales, which can be inferred from 
the study of disc frequencie s (e. g.lHaisch et al.l (120011 ) ). are of order a few million years 
(e. g. lHernandez et all fl2007t ): iMamaiekl ( 12009( 1 ). Dispersal t ime-scales, as determined from 
the study of infrared colours of young stars, appear to be much shorter, indicating that the 
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dispersal mechanism m ust be fast (e. g. lKenyon fc Hartmannl fl 19951 ): iLuhman et al.l ( 120101 ) : 



Ercolano et ID fl2011bh ^. Such observed two-timescale behaviour has favoured the develop 



ment of disc dispersal models that involve a rapid disc clearing phase, contra ry to the pre- 



dictions of simple viscous drai nin g, and in agreement with photoevapor a tion dClarke et al. 



2001; Alexander et al.. 2006a b; Ercolano. 2008; Ercolano et al.. 2009; Gorti et al.. 2009 



Owen et al.l . l2010l . l2011aU bl. 120121 ) or possibly planet formation ( Armitage fc Hansen! Il999l ). 



The interpretation of infrared colours in relation to the evolutionary state of a disc is, how- 
ever, far from being trivial. This is particularly true with regards to the classification of 
transition discs, the latter being intended as objects caught in the act of disc dispersal. The 
evolution of the dust component in a disc is mirrored by the evolution of colours in the 
infrared pl ane. By means of radiative transfer modelling, Ercolano. Clarke. Sz Hjdj (2011, 
henceforth lECHllI ) identified the regions in the K - [8] vs. K - [24] plane where primordial 
discs, discs with inner-holes (i. e. presumably being dispersed from the inside-out) and discs 
which lose mass homogeneously at all radii, are expected to be found. Their study, which 
was limited to M-stars, showed that in the case of the cluster IC348, most discs disperse 
from the inside-out and undergo the transition on a short time-scale, as predicted by stan- 
dard phot oionisat i on m odels. These conclusions are in contrast with the conclusions of 



Currie fc Kenyonl (120091 ). who claimed instead a large number of 'homogeneously depleting' 



discs, for the same cluster. Such discrepancies highlight the need for detailed modelling in 
the interpretation of IR colours of discs. The study of lECHllI was restricted to M-stars 
in only one cluster, which prevented the authors from being able to make a more general 
statement with regard s to disc dispersal. Here we present results fro m a fort hcoming paper 



(IKoepferl et al.l . 120121 ) . which significantly improves on the work of lECHllI by performing 



further RT calculation to evaluate evolutionary tracks in the IR colour plane for stars of 
different spectral types. We then apply our results to the photometric data of 15 nearby 
star-forming regions, that we collected from the literature, in order to address the question 
of what is the preferred mode of disc dispersal. 



2 Disc evolution in the Infrared two colour plane 

By means of radiative transfer modelling we identify 5 areas in the K - [8] vs. K - [24] plane 
corresponding to discs at different evolutionary stages or with different geometries, these are 
shown in Figure [TJ Namely, primorial optically thick (flared and/ or mixed) discs occupy 
area A, while disc-less objects cluster in area B and area C is for primorial ultra-settled (flat) 
discs. Discs that are in the act of dispersal belong to area D and E, where the former is for 
discs that are clearing from the inside-out, and the latter is for discs that are progressively 
going optically thin homogeneously at all discs radii (as would be expected from viscous 
evolution alone). Our classification scheme allows us to use currently available infrared pho- 
tometric surveys from nearby star-forming regions in order to address the question of what 
is the preferential mode of disc dispersal. We therefore applied our colour-colour diagnostic 
diagrams to classify 1529 objects in 15 nearby star-forming r egions. As an example, we 



overplot the data for Taurus (taken from ILuhman et al.l (120101 )) to the diagnostic diagram 



in Figure [TJ In summary, 39 % of the objects out of the entire sample lie in the primor- 
dial disc region whereas between 31 % and 32 % disperse their discs from the inside-out and 
up to 22 % of the sources have already lost their disc. So, almost a third of the available 
sources are currently clearing their discs from the inside-out. Less than 2% of the objects 
lie in the homogeneous draining region E. We interpret this result as strong evidence against 
homogeneous disc depletion as the main disc dispersal mode. 
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Taurus (Luhman) 



G0<SpT<K4.5 K4.5<SpT<M2.5 




K-[24](mag) K-[24](ma9) 



Figure 1: Disc evolution diagnostic dia- 
gram applied to the YSOs in the 1 Myr 
old cluster Taurus. 20 sample points lie 
outside the limits of this plot, but they are 
still included in the final statistics. (51% 
primordial optically thick (A), 0% pri- 
mordial ultra-settled (C), 23% disc- less 
(B), 14% inside-out clearing (D), < 1% 
homogeneous draining (E)). 



3 Dispersal time-scales across spectral types 



With the sample of YSOs in different star-forming regions becoming larger and larger, it is 
now possible to estimate the typical time-scales for the disc dispersal phase, even though 
cluster ages of course always introduce a large uncertainty in the estimates. The disc- 
evolution time-scale of a star-forming region can be roughly estimated by multiplying the 
age of the region by the ratio of the number of evolved object to the total number of objects 
in the region. Typical transition time-scales for the considered star-forming regions are of 
order 10 5 yrs. The average time-scale across all spectral types is 6.9 • 10 5 yrs and roughly 
the same as the average cluster time-scale of 6.6 ■ 10 5 yrs. This is partially because there 
is no significant difference for timescales amongst spectral types. We further illustrate that 
disc dispersal time-scales appear to be independent of spectral type. We plot in Figure [21 
the time-scale ratio between K and M-stars and show that this ratio is consistent with unity. 
This suggests that there is no si gnificant dependence o f the time-scale on stellar mass, as 
has already been pointed out by lErcolano et al.l ( 2011al ). who performed a spatial analysis 
of the distribution of K and M-stars with discs in young star-forming regions and found no 
significant difference in the distributions. 



IC34B Taurus Taurus Upper Sco OBIbl 25 Ori NGC1333 NGC2264 Cha I 
(Luhman| (Rebull) 

star-forming regions 



Figure 2: Time-scale ratio r(ii)/r(iii) for K 
and M stars for the star-forming regions 
considered in this work with N to t > 10. 
Three star-froming regions are not listed, 
because they lack evolving objects in one 
or both spectral type intervals. 
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4 Summary 

We have calculated the SEDs of protoplanetary discs of different spectral types, geometries, 
settling and inclination. We then considered the evolution of the infrared colours (K - [8] 
vs. K- [24]) of the model disks as they disperse according to different scenarios (homologous 
depletion, inside-out and outside-in clearing). Based on our models we propose a new diag- 
nostic infrared colour-colour diagram to classify the evolutionary stage of YSOs. We have 
applied our infrared colour-colour diagnostic diagram to classify YSOs in 15 nearby star- 
forming regions and study the evolution of their disc populations. We estimate time-scales 
for transition phase of typically a few 10 5 years independent of stellar mass. We conclude 
that, regardless of spectral type, current observations point to fast inside-out clearing as the 
preferred mode of disc dispersal. 
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